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NOMENCLATURE 


CI  ' n.  u An 

1 00  U 


-M 


collecting  disc  area 


ion  flux  coefficient 


diameter 
electric  field 
charge  on  electron 
mobility 

Boltzmann's  constant 


reduced  mass 
number  density 
pressure 

collision  cross  section 


speed  ratio 


temperature 

velocity 

thermal  velocity 
drift  velocity 


r0/re 

+</*e 


P 


mean  free  path 
Debye  length 
polarizability 


v 


a 


particle  diameter 
electric  potential 


Subscripts 


free  stream  conditions 
electron 
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A unique  computational  procedure  to  predict  the  flux  of  positive  ions  to 
a rocket-borne  mass  spectrometer  has  been  developed  under  previous  contracts. 
This  procedure  is  the  Direct  Simulation  Monte  Carlo  Method  applied  to  charged 
particles.  The  method  developed  initially  for  neutral  gases  has  been  extended 
to  a plasma  in  the  presence  of  a charged  botfy  with  a consistent  electric  field. 
Since  the  electric  potential  distribution  required  a solution  of  Poisson's 
equation,  the  charged  particle  calculations  required  an  iteration  procedure 
between  the  Monte  Carlo  results  and  a finite  difference  solution  of  Poisson's 
equation.  The  detailed  formulation  and  results  obtained  for  the  collection 
of  positive  ions  are  given  in  references  1 through  4. 

In  the  present  study,  the  computer  code  developed  for  the  positive  ion 
collection  has  been  modified  to  predict  the  flux  of  negative  ions  to  a rocket- 
borne  mass  spectrometer.  The  major  modification  was  the  calculation  of  the 
electron  trajectories  which  were  ignored  in  the  previous  studies. 


* 
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2.  TECHNICAL  APPROACH 

2.1  PROBLEM  FORMULATION  FOR  NEGATIVE  ION  ANALYSIS 

In  previous  studies,  e.g.,  references  1 through  3,  the  TRW  Monte  Carlo  direct 
simulation  has  been  used  to  predict  the  positive  Ion  flux  to  a rocket  borne  mass 
spectrometer.  A systematic  variation  of  parameters  was  performed  for  the  follow- 
ing variables:  speed  ratio,  mean  free  path,  Debye  length,  and  applied  collecting 
potential.  The  predictions  and  trends  resulting  from  the  parametric  study  will 
be  helpful  in  understanding  flight  data  and  in  the  design  of  future  instruments. 
Also  included  in  the  previous  studies  were  the  following  effects:  internal 
degrees  of  freedom,  angle  of  attack,  and  space  charge  build-up  on  the  vehicle 
sidewalls. 

The  assumptions  and  the  complete  formulation  of  the  Monte  Carlo  simulation 
are  detailed  in  Reference  1.  Most  important  of  the  assumptions  regards  the 
treatment  of  the  electrons.  For  the  positive  ion  collection,  the  front  face  is 
negatively  charged  which  repels  the  electrons  and  attracts  the  positive  ions. 

The  negative  front  face  potential  was  assumed  large  enough  so  that  the  electrons 
could  be  taken  to  be  in  equilibrium  with  the  local  potential  and  determined  by 
the  Boltzmann  distribution. 


Ilf-  = e(e*/kT  ) 

FT  e e 

e« 

where 

Ng  = electron  number  density 

Na  = free  stream  value 
e°° 

♦ = electric  potential 

Te  = electron  temperature 
e = charge  on  electron 

This  assumption  was  shown  to  be  accurate  in  Reference  1 and  was  then  used  through- 
out the  previous  studies.  The  great  simplification  of  this  assumption  was  that 
the  electron  trajectories  could  be  ignored  and  only  the  ion  motion  was  computed 
with  the  only  coupling  between  the  ion  and  electron  distributions  occurring 
through  the  electric  potential  determined  from  Poisson's  equation 
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where  * = * non-dimensional  electric  potential 

oo 

x = mean  free  path 

ao 

x = Debye  length 

U -f 

+ V 

^ * jj—  * positive  Ion  number  density  ratio 

i GO 

- Ne 

Ne  = jj—  * electron  number  density  ratio 

poo 

In  the  present  study,  the  negative  ion  flux  to  a rocket  borne  mass  spec- 
trometer is  investigated.  The  collection  of  negative  ions  requires  a positive 
collecting  potential  on  the  front  face.  This  positive  potential,  however,  also 
attracts  the  electrons  which  will  reduce  the  collecting  potential  and  the  effec- 
tiveness of  the  instrument.  The  sampling  of  negative  ions  requires  a circular 
double  disk  configuration;  Reference  5,  on  the  front  face  to  repel  the  electrons 
while  still  collecting  negative  ions.  The  geometry  and  the  potential  distribution 
along  the  front  face  of  the  instrument  payload  will  be  approximated  as  shown  in 
Figure  1.  Although  the  face  potential,  4>o,  would  be  greater  than  the  mask  poten- 
tial, 4>e,  the  mask  has  greater  area.  Therefore,  only  the  proper  ratio  of  4>o  to 
4>e  will  allow  the  ions  to  pass  through  the  collecting  orifice  for  a given  set  of 
free  stream  and  boundary  conditions.  The  electron  motion  in  this  case  cannot  be 
ignored,  and  the  equilibrium  assumption  for  the  electrons  is  no  longer  valid. 

This  is  the  major  complication  introduced  by  the  negative  ion  collection  which 
was  not  encountered  in  the  sampling  of  positive  ions. 

These  complications  can  be  seen  by  comparing  the  electron  and  ion  velocities. 
The  thermal  velocities  for  the  electrons  and  ions  are  given  by 


v,  = y1 2 ^Ti  for  ions 
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for  electrons 
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♦o  = positive  collection  potential  for  negative  ion 
$e  = negative  electron  mask  potential 

Figure  1.  Front  Face  Geometry 

The  ratio  of  these  speeds  are  given  by 


which  for  equal  temperatures,  = Tg , and  an  ion  molecular  weight  of  30  becomes 

!l=235 


Since  the  vehicle  or  free  stream  velocity,  Uw , is  roughly  3 to  5 times  the  ion 
thermal  speed,  the  ions  are  characterized  by  the  free  stream  velocity 


Ui  a U-  “ 4 ?1 

whereas  the  typical  electron  velocity  is  the  thermal  velocity 


ue  ^ Ve  a 235  a 60 


Therefore,  for  a given  time  increment,  the  electron  travels  a distance  60  times 
the  ion  displacement.  If  the  time  step  in  the  Monte  Carlo  calculations  is  given 
by  At . for  the  ions,  the  result  requires  60  electron  time  steps  to  equal  At^. 

In  other  words,  each  ion  calculation  requires  60  electron  calculations.  1 hi s 
number  can  be  reduced  by  recalling  that  both  the  electrons  and  ions  are  trace 
species  and  that  encounters  between  charged  particles  are  neglected  except  through 
their  coupling  in  the  electric  potential.  The  only  direct  interactions  are  col- 
lisions between  the  charged  particles  and  the  neutral  molecules  which  can  be 
characterized  by  the  mean  free  paths 
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for  ions/neutrals 


where  Nn  = number  density  of  neutrals 

Qen»  Qin  = collision  cross  sections 

Assuming  hard  spheres,  the  collision  cross  sections  for  two  different  particles 
is  given  by 


where 


al » °2  = ParT1cle  diameters 
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For  the  ions,  & an  and  Q^n  a it  , but  for  the  electrons  og  « on  and 

Pena  * an  • The  ratio  mean  free  paths  is  then  given  by 


Therefore,  if  an  ion  travels  a distance  in  time  increment  At  an  electron 
displacement  can  be  xgn  ^ 4 x^n.  This  reduces  the  number  of  electron  time  steps 
to  equal  At^  from  60  to  60/4  = 15.  In  rough  terms,  this  means  that  the  negative 
ion  collection  calculations  will  take  approximately  15  times  as  long  as  the 
equivalent  positive  ion  calculation  which  ignored  the  electron  trajectories. 

The  procedure  to  be  used  in  the  negative  ion  sampling  and  a number  of  ways 
to  reduce  the  computing  time  is  discussed  In  the  next  section. 

2.2  COMPUTATIONAL  PROCEDURES 

The  calculation  procedure  for  the  collection  of  negative  ions  is  shown  in 
the  flowchart  on  Figure  2.  The  calculation  is  split  into  two  parts  to  reduce 
the  computer  time  requirements  by  eliminating  unnecessary  repetition  and  proceeds 
in  the  following  steps: 

(1)  Initialize  the  problem  by  specifying  the  geometry  and  flow  condi- 
tions, e.g. , 

• Vehicle  velocity 

• Mean  free  path  for  neutral  molecules 

• Debye  length 

• Electric  potential  distribution  on  the  vehicle. 

(2)  Since  the  ions  and  electrons  are  trace  species,  the  neutral  molecule 
flowfield  can  be  determined  without  considering  the  ions  and  electrons. 
The  neutral  flowfield  computed  by  the  Monte  Carlo  direct  simulation 
method  Is  used  as  a scattering  field  for  the  ion  and  electron  distri- 
butions. This  calculation  Is  performed  only  once  for  a given  mean 
free  path  and  can  be  used  for  all  Debye  lengths.  In  the  collisionless 
limit,  this  calculation  is  bypassed. 

(3)  The  electric  potential  is,  in  general,  determined  from  the  Poisson 
equation  given  by 


'•  • fe)  <V  - V - V 


where  all  the  variables  are  identical  to  those  defined  in  Equation  (2) 
with  the  addition  of 

- - Ni" 

N.  = = negative  ion  number  density  ratio 

^ 1«0 

At  this  point,  the  ion  and  electron  number  densities  have  not  been 
determined.  Therefore,  the  first  interation  for  the  electric  potential 
is  found  by  solving  the  Laplace  equation 


r*  = 0 


which  would  be  exact  if  the  Debye  length  were  infinite,  e.g.,  x^  = °°. 

In  the  general  case,  the  Laplace  solution  serves  as  the  first  iteration. 

(4)  4a  and  4b  --  The  ion  and  electron  flowfields  are  computed  by  the  Monte 
Carlo  method  for  the  neutral  field  determined  in  step  (2)  and  the 
potential  distribution  found  in  step  (3).  Note  that  these  flowfields 
provide  the  solution  for  a given  mean  free  path  and  a Debye  length  of 
infinity.  It  has  been  found,  however,  in  previous  calculations  for 
positive  ions  that  the  Laplace  solution  to  Poisson's  equation  occurs 
for  Debye  numbers,  xD^D,  of  order  unity.  Therefore,  these  solutions 
obtained  at  this  step  are  probably  applicable  for  Debye  numbers, 

xD^0  >1.  It  has  also  been  shown  in  Reference  4 that  a significant 
region  of  the  ionosphere  of  interest  for  ion  sampling  is  characterized 

by  XD/D  > 1 * 

(5)  After  the  individual  number  densities  are  determined,  the  charge 
density 

p = 1T.|+  - TTj“  - TTe 

which  appears  on  the  right  hand  side  of  Equat1on(3)  is  evaluated  through- 
out the  flowfleld. 

The  calculations  from  (1)  through  (5)  are  not  repeated  In  the  calculation 
for  a given  set  of  flow  parameters  and  are  separated  for  subsequent  calculations 
and  designated  as  Part  I.  This  division  is  also  made  since  steps  (1)  through  (5) 
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I 

H 

are  performed  for  a Debye  length  of  = « which  must  be  repeated  in  principle 
for  every  calculation  where  all  the  parameters  are  constant  except  the  Debye 
length.  Therefore,  to  eliminate  repetition  of  a time  consuming  procedure,  the 
information  generated  in  Part  I,  e.g., 

(1)  Boundary  and  initial  conditions 

(2)  Neutral  molecule  flowfield 

(3)  Laplace  solution 
(4a)  Ion  flowfield 
(4b)  Electron  flowfield 

(5)  Charge  density 

will  be  stored  on  a parmanent  file  and  recalled  and  used  for  calculations  involv- 
ing different  Debye  lengths  for  the  same  values  of  the  other  parameters. 

Part  I I. a which  is  dependent  on  the  Debye  length  is  shown  in  Figure  3 and 
begins  with  step  (6). 

(6)  A solution  to  Poisson's  equation  for  the  charge  density  distribution 

j 

computed  in  step  (5)  and  the  proper  Debye  length,  xQ. 

, The  following  steps  describe  the  basic  iteration  procedure  to  determine  the 

electron  and  ion  flowfields  consistent  with  the  Poisson  equation  and  the  neutral 
molecule  distribution  from  step  (1).  Since  this  is  a new  computation  procedure, 
the  optimum  method  to  perform  these  calculations  has  not  been  determined.  There- 
fore, two  methods.  Part  I I. a and  Part  Il.b  were  developed  and  compared.  Note 
that  in  either  case  the  ion  and  electron  motions  are  coupled  only  through  the 
electric  potential  since  both  are  trace  species  and  only  encounters  with  the 
neutral  background  molecules  are  important.  The  first  method  is  shown  in 
Figure  3 in  steps  (7)  through  (10). 

(7)  Electron  movement  — The  ions  and  neutrals  are  assumed  stationary 

in  this  step  and  the  electron  motion  is  computed  for  the  current 

electric  potential  distribution.  As  discussed  in  the  previous 

section,  the  higher  speeds  associated  with  the  electrons  requires 

that  approximately  60  time  steps  are  required  for  the  electrons 

for  each  ion  time  step.  However,  the  difference  in  the  man  free 

paths  reduces  this  ratio  from  60  to  15.  Therefore,  for  a fixed 

movement  time  Interval,  At  , the  electrons  require  15  more  cai- 
rn 

culations  than  for  the  ions. 
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POISSON  SOLUTION 
FOR  THE  PROPER  OEBYE  LENGTH,  An 
INITIAL  SOLUTION  IN  ITERATION 


Figure  3.  Flow  Chart  for  Part  II. a;  Iteration  Procedure  for 
Negative  Ion  Collection 
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(8)  Ion  movement  — The  electrons  and  neutrals  are  kept  stationary 
and  the  ion  motion,  both  positive  and  negative,  is  computed  for 
the  current  electric  potential.  Note  that  the  electric  potential 
has  not  been  updated  between  (7)  and  (8),  but  the  influence  of 
recalculating  the  potential  after  step  (7)  will  be  investigated. 

(9)  Electric  potential  — The  new  charge  density  and  electric  potential 
will  be  computed  for  the  electron  and  ion  distribution  calculated 
in  steps  (7)  and  (8). 

(10)  Convergence  test  — Convergence  of  the  electron  and  ion  distributions 
for  a consistent  Poisson  solution  and  the  neutral  molecule  distribution 
will  be  assessed  in  step  (10).  This  will  be  accomplished  by  checking 
the  charged  particle  distributions  for  the  current  results  with  the 
previous  iteration.  If  the  comparisons  are  satisfactory,  the  results 
are  printed;  if  not,  the  iteration  continues  by  returning  to  step  (7). 

The  alternate  procedure  to  determine  the  electron  and  ion  distributions,  labeled 
Part  II. b,  is  shown  in  Figure  4.  In  both  Part  II. a and  II. b,  the  initial  calcu- 
lations performed  in  steps  (1)  through  (6)  are  identical  and  only  the  iteration 
procedure  is  changed.  Part  II. b has  three  iteration  loops  compared  to  one  for 
Part  1 1. a and  is  described  below. 

(7)  Electron  movement  and  iteration  --  The  ion  and  neutral  distributions 
are  fixed  and  the  electron  flowfield  is  determined.  This  calculation 
differs  from  Part  II. a in  that  the  "steady  state"  electron  distribution 
is  determined  for  a Poisson  solution  consistent  with  the  electron 
motion  and  the  stationary  ion  and  neutral  fields.  Since  only  the 
electrons  are  in  motion,  the  time  steps  are  independent  of  the  ion 
time  steps  and  are  not  constrained  as  in  step  (7)  of  Part  II. a. 

(8)  Ion  movement  and  iteration  — Using  the  electron  density  and  the 
potential  found  in  step  (7),  the  positive  and  negative  ion  trajectories 
are  determined.  Again,  the  neutral  molecules  are  fixed  and  the 
Poisson  solution,  consistent  with  the  ion  motion,  is  determined. 

(9)  Convergence  test  — The  electric  potential  distributions  computed 
in  (7)  and  (8)  are  compared  for  convergence.  If  the  potentials  are 
acceptably  similar,  the  results  are  printed;  if  not,  the  iteration 
returns  to  step  (7). 

At  the  start  of  the  present  study  no  direct  comparisons  were  available  between 
the  procedures  and  potential  advantages  are  apparent  in  both.  For  example,  in 
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Figure  4.  Flow  Chart  for  Part  Il.b;  Iteration  Procedure 
for  Negative  Ion  Collection 


L.t. 


Il.b  there  is  no  constraint  on  the  time  step  for  the  movement  of  the  electrons 
as  in  II. a.  Therefore,  the  electron  iterations  should  be  as  fast  as  the  ion 
iterations,  but  the  convergence  between  the  ion  and  electron  fields  could  elim- 
inate any  potential  advantage  over  II. a.  In  the  method  proposed  in  II. a,  the 
electron  iteration  was  estimated  to  be  15  times  greater  than  the  ion  calculation 
for  a given  time  interval.  Therefore,  this  method  is  not  feasible  unless  the 
electron  time  step  can  be  increased  in  practice.  For  example,  if  a typical  ion 
convergence  without  electrons  takes  10  minutes  on  the  CDC  6600,  the  electron 
calculations  would  take  150  minutes  which  is  clearly  unacceptable  (for  low 
altitude  results  involving  many  collisions,  this  number  could  be  doubled).  An 
upper  limit  on  the  entire  calculation  of  approximately  60  minutes  requires  that 
the  time  factor  be  reduced  from  15  to  5.  Another  computational  method  which 
would  accelerate  the  calculations  would  be  to  make  better  use  of  weighting  factors 
for  the  charged  particles  so  that  the  number  of  particles  in  the  calculation  could 
be  reduced.  No  attempt  to  implement  these  efficiencies  was  made  in  the  previous 
calculations  since  the  computation  times  for  the  positive  ion  collection  were 
well  within  reasonable  constraints.  Both  methods  have  been  implemented  and  a 
few  test  cases  have  been  computed.  From  these  preliminary  results  the  method 
designated  as  Il.b  appears  to  converge  faster.  No  conclusions  were  made  during 
the  reporting  period  due  to  the  observation  that  the  electron  effects  were  not 
simulated  properly.  This  is  discussed  in  the  next  section. 

2.3  CHARGED  PARTICLE-NEUTRAL  INTERACTION 

Monte  Carlo  calculations  performed  for  the  electron  motion  in  a neutral  gas 
at  densities  where  collisions  are  Important  indicated  that  the  drift  velocities 
were  too  high.  After  confirming  the  correctness  of  the  numerical  procedure,  the 
basic  assumptions  for  the  simulation  were  reexamined. 

In  the  past  simple  hard  spheres  were  used  for  the  charged  particle-neutral 
scattering  model  and  the  collision  cross  section  was  computed  from 


where  ap  * charged  particle  diameter 
on  • neutral  particle  diameter 
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For  the  ions  (p 


end  for  the  electrons  (p  a e) 


The  simple  hand  sphere  model  used  for  convenience  in  the  earlier  studies  was  justi- 
fied on  the  basis  that  the  gross  features  of  the  ion  collection  process  were  not 
clearly  understood  and  in  the  collection  of  positive  ions  the  electron  trajectories 
were  not  computed.  However,  the  present  task  of  simulating  the  collection  of  neg- 
ative ions  requires  accurate  calculations  of  the  electron  motion. 

The  proper  description  of  the  charged  particle-neutral  encounter  can  be  rep- 
resented by  an  Interaction  potential  of  the  form 

V (r)  * Ae“*r  + Cnr*n 
n * 4 


where  the  first  term  represents  strong,  short  range  repulsive  forces  due  to 
ordinary  neutral  particle  gas  kinetics  and  the  succeeding  terms  represent  attrac- 
tive forces  which  arise  between  the  charged  particle  and  the  polarized  structure 
of  the  neutral. 

In  the  calculations  presented  here  a 4th  power  attraction  suggested  origi- 
nally by  Langevin  was  used.  The  attractive  force,  due  to  polarization  of  the 
neutral  molecule  is  given  by 


where  u = polarizability 

e = electronic  charge 
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The  current  Monte  Carlo  calculations  treat  the  charged  particle-neutral  encounter 
as  a binary  collision  involving  both  attractive  and  repulsive  forces.  The  particle 
trajectories  due  to  the  repulsive  force  can  be  separated  into  two  effects:  (1) 
Isotropic  scattering  at  an  impact  parameter  inversely  proportional  to  the  relative 
velocity  of  the  collision,  (2)  isotropic  scattering  by  a hard  sphere  core  of  fixed 
radius.  The  dominance  of  (1),  (2)  or  the  attractive  force  depends  on  the  value 
of  the  Impact  parameter.  For  low  energy  encounters  the  attractive  force  would 
dominate,  whereas  for  high  energy  encounters  strong  repulsive  forces  would  be 
most  Important  (see  Figure  5). 

The  collision  diameter  for  the  attractive  force  is  found  to  be 

o 


where  m = 


m„  + in 
p n 


uf  ■ relative  velocity  of  the  encounter 


This  diameter  can  be  compared  to  the  standard  hard  sphere  diameter  on  to  illustrate 
the  polarization  effect.  For  example,  oxygen  has  a polarizability 

v - 1.6  x 10"24  cm3 


and  a hard  sphere  diameter 


on  » 3 x 10“®  cm 


The  reduced  mass  for  the  ion-neutral  interaction  becomes 


* « m_ 

_ _ n 1 'x*  n 


ftij' +_mn"  - r (m1  ' V 


At  300° K the  relative  velocity  can  be  approximated  by 


2 * 2kT 


1.4  x 10 


9 cnr 
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Finally  the  ratio  of  the  diameters  is  found  to  be 


This  result  implies  that  the  collision  cross  section  for  oxygen  ion-neutral  encoun- 
ters at  300°  K is  approximately  25  times  the  typical  hard  sphere  value. 

To  investigate  the  results  of  this  collision  model  a simple  relaxation 
process  has  been  studied.  Consider  a weakly  ionized  gas  composed  of  electrons, 
positive  ions  and  neutral  particles  which  are  initially  in  equilibrium  at  a single 
temperature.  The  positive  ions  are  singly  charged  and  of  mass  equal  to  the  neutral 
particles.  A prescribed  constant  electric  field  is  instantaneously  applied  and 
subsequently  maintained  and  the  molecular  motion  of  the  gas  is  computed  in  time 
as  a steady  state  is  approached.  The  physical  situation  to  which  this  corresponds 
could  be  Imagined  to  consist  of  a gas  contained  between  parallel  electrodes  of 
unlimited  extent  and  great  separation.  Attention  is  focused  upon  a thin  slab  of 
this  gas  within  which  the  spatial  gradients  can  be  ignored.  Since  all  charged- 
charged  interactions  are  contained  in  the  electric  field,  which  is  known  and 
constant  for  this  case,  the  electron  and  ion  motions  are  independent.  Only  the 
encounters  of  the  charged  particles  with  the  neutral  particles  while  under  the 
influence  of  the  electric  field  determine  the  behavior  of  the  charged  particles. 

Quantities  of  interest  for  this  study  are  the  following: 

1)  the  steady  state  drift  velocity  per  unit  field,  e.g.  the  mobility 

2)  the  relaxation  times 

3)  the  temperature  and  shape  of  the  distribution  function  in  the  steady 
state 

In  the  ion  collection  process  by  a rocket  borne  mass  spectrometer  biased  at  a 
potential  $ the  parameter  e$/kT  is  used  to  characterize  the  strength  of  the 
electric  field.  In  the  present  example  of  a gas  contained  between  parallel 
electrodes  the  parameter  E/P,  field  strength  over  neutral  gas  pressure,  is  used. 

The  acceleration  of  a particle  with  charge  e and  mass  m in  an  electric  field  E is 
eE/m.  The  energy  acquired  by  the  particle  from  the  field  between  collisions  is 
approximately  eE*x  where  x is  the  mean  free  path.  Since  x is  inversely  propor- 
tional to  the  density 
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E-A  ^ E* 


E 

F 


T 


or 


E 

P 


Therefore  the  ratio  of  the  field  energy  to  thermal  energy  is  proportional  to  E/P. 
This  result  assumes  the  field  strength  is  low  and  that  the  masses  of  the  ion 
and  neutral  particle  do  not  differ  greatly  (thus  excluding  electrons).  For 
arbitrary  field  strength  and  ion  mass,  E/P  is  less  useful.  However,  most 
experimental  data  are  presented  in  terms  of  it. 

Results  of  computations  made  with  this  model  is  shown  in  Figures  6 through 
9.  Here  the  product  of  the  mobility 


K = 


where  = drift  velocity 
E = electric  field 

and  the  pressure,  p,  in  atm  is  shown  for  the  electrons  in  a neutral  gas  with 
molecular  weight  29.  The  calculations  compare  adequately  with  the  experiments 
which  have  been  influenced  by  effects  not  modeled  in  the  simulation  such  as 
charge  transfer,  quantum  mechanical  resonance  effects  and  clustering. 

2.4  NEGATIVE  ION  COLLECTION  - COLLISIONLESS  RESULTS 

A series  of  calculations  have  been  performed  to  predict  the  collection  of 

negative  ions  to  the  probe  geometry  shown  in  Figure  1.  The  ratio  of  the  positive 

collection  potential  to  the  negative  electron  mask  potential,  y = *0/*e,  was 

fixed  at  -10  and  the  ratio  of  the  disk  radii,  a = r /r  , was  varied.  The 

(c\  o e 

critical  value  of  y is  given  by'  ' 


For  yalues  of  |y|<|yc  a repelling  potential  exists  along  the  centerline  which 
retards  the  collection  of  the  negative  ions.  This  function  is  shown  in  Figure  10 
and  Indicates  that  for  y = -10  a repelling  potential  will  exist  for  a < 0.4. 
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Figure  5.  Charged-Neutral  Scattering  Model 


Figure  6.  Electron  Mobility 


ure  7.  Electron  Velocity  Distribution 


SAMPLE  SIZE 


Figure  11  shows  the  variations  of  the  centerline  potential  for  the  following 
conditions, 

Speed  ratio  s = 3.10  = 200  *e  = -20 

Xoo 

fe;  sen  number  = 1000 


XD 

Debye  number  = 


1000 


and  a from  0.08  to  0.75.  The  potentials  for  a = 0.08  and  0.25  exhibit  the  change 
in  sign  of  potential  as  expected  and  the  result  for  a - 0.42  appears  to  be  at  the 
transition  point.  For  the  largest  value  of  a = 0.75  the  potential  decays  in  a 
monotonic  way.  The  flux  of  positive  ions  collected  by  the  probe  point  is  shown 
in  Figure  12  as  a function  of  the  collecting  radius.  For  a = 0.08  which  resulted 
in  a maximum  repelling  potential  of  approximately  <j>  = -7  there  is  no  flux  of 
negative  ions  to  the  front  face.  However,  for  a = 0.25  with  a small  repelling 
potential  of  $ = -1  the  kinetic  energy  of  the  negative  ions  was  sufficiently 
higher  than  the  retarding  electric  field  to  establish  a flux  coefficient  at  the 
stagnation  point  of  approximately  Cj  = 10.  As  the  value  of  a is  increased,  the 
flux  coefficient  changes  only  slightly  as  a varies  from  0.58  to  0.79. 

Although  these  calculations  were  performed  for  a Debye  number  large  enough 
to  neglect  the  charge  density  effects,  the  electron  flow  field  and  flux  to  the 
front  face  must  be  computed  to  determine  the  current. 
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3.  CONCLUSIONS 


The  efforts  to  extend  the  Monte  Carlo  Method  to  the  collection  of  negative 
ions  have  been  shown  to  be  more  difficult  than  anticipated.  It  was  found  that 
the  electron  trajectories  must  be  determined  carefully  and  that  the  effect  of 
polarization  are  important  in  electron-neutral  collisions.  In  addition,  the 
Monte  Carlo  computer  code  which  was  developed  for  the  negative  Ion  collection 
stretched  the  computer  resources  of  the  CDC6600.  Future  efforts  will  include 
the  implementation  of  more  efficient  programming  techniques  to  reduce  the  core 
requirements  and  computation  times. 
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